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We compare Monte Carlo simulations of knotted and unknotted polymers whose ends are con- 
nected to two parallel walls. The force / exerted on the polymer is measured as a function of 
the separation R between the walls. For unknotted polymers of several monomer numbers N, the 
product fN" is a simple function of R/N v , where v ~ 0.59. By contrast, knotted polymers exhibit 
strong finite size effects which can be interpreted in terms of a new length scale related to the size 
of the knot. Based on this interpretation, we conclude that the number of monomers forming the 
knot scales as AT*, with t = 0.4 ± 0.1. 
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Entanglements are unavoidable in long polymers and 
influence their properties jj] . Knots are found in proteins 
, and present an obstacle that needs to be overcome 
in the transcription of DNA S . An increasing number 
of experiments can now probe the detailed properties of 
knotted molecules Micro-manipulation techniques 
H enable direct measurements of mechanical properties 
of a single molecule, and it is even possible to probe the 
behavior of artificially knotted DNA ||. However, in- 
corporation of topological constraints into the statistical 
mechanics of polymers |?J remains a difficult theoretical 
challenge since the resulting partition of phase space into 
accessible and inaccessible regions cannot be easily im- 
plemented. Nevertheless, some progress has been made 
in understanding the role of knots in loop polymers; e.g. 
the relative probabilities for appearance of different knots 
in self-avoiding (SA) loops || has been characterized. 
However, much less is known about the typical shapes 
and physical properties of knots. 
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FIG. 1: Two configurations of a knotted chain of N — 335 
bonds, with its ends attached to parallel walls separated by 
50a. The walls are not shown and the chains were rotated 
for clarity of view. The (grey) knotted regions consist of 32 
bonds in (a) and 112 bonds in (b). 

As a topological feature, a knot can be rigorously iden- 
tified only by specifying the entire shape of a closed chain 
(ring). However, it is natural to identify a segment where 
the knot is located, and consequently to talk about its 
size and statistics. It is natural to pull on the ends 
of a string to see if a small knot remains in the mid- 
dle. Our eyes tend to identify knotted segments in this 
manner, as exemplified by the grey-shaded monomers of 
the polymers in Fig. [I]. While pulling on a string makes 
tight knots, the question of whether in a random unforced 
chain the knot is spread over the whole curve, or localized 
(tight) on a small portion, is still not fully resolved. Sev- 
eral recent works |j show that knots in SA walks confined 
to a two-dimensional (2D) plane are strongly localized, 
i.e. the mean number of monomers in the knotted region 
Nk does not depend on the total number of monomers 
N. Studies of three-dimensional (3D) knotted polymers 
are hampered by the difficulty of identifying the knotted 
region. Katritch et al. [[l0| examined the size distribu- 
tion of knots in 3D random rings, by removing segments 
of the ring, attaching them to infinite straight lines, and 
checking if the resulting structures were knotted. Al- 
though the method has a certain probability to fail (i.e., 
the procedure itself may create or remove a knot), it nev- 
ertheless suggests that knots in such rings are localized. 
By contrast, some studies jnj had earlier indicated that 
for moderate N, the radius of gyration R g of a knot is 
strongly influenced by its complexity, leading to the con- 
clusion that the knots might be spread out over the entire 
loop p2[ . More recent numerical results [^3[ |l4| provide 
evidence that R g is asymptotically independent of the 
knot type, hinting that knots are localized. 

Since correlation functions of SA walks are power laws, 
it is reasonable to expect that the distribution of the 



2 



number of monomers n forming a knot is also a power 
law, i.e. p(n) — Cn l ~ 2 . For t < 1 the normalization 
coefficient C is determined by the microscopic (short dis- 
tance) properties of the chain, while for < t < 1, the 
expectation value of the mean knot size Nk depends on 
the total length, growing as N t . In such a case, we say 
that the knot is weakly localized. For t < 0, the knot 
is strongly localized in the sense that Nk is determined 
by the microscopic cutoff. In this paper, we attempt 
to quantify the tightness of knots in 3D polymers by 
comparing the force-extension relations of knotted and 
unknotted chains. As in Ref. [jl4j, we find similarities 
between the statistical properties of knotted chains and 
those of shorter unknotted chains, a consequence of the 
fact that knotted segments are statistically denser than 
unknotted ones. Using the reduction in the effective num- 
ber of monomers as an operational definition of knot size, 
we find N k ~ N*, with t = 0.4 ± 0.1. 

Much of the current understanding of the scaling prop- 
erties of long polymers is based on renormalization group 
ideas |l5). The N — > oo limit in polymers corresponds 
to approaching a fixed point similar to that describing 
criticality in thermal phase transitions. In particular, it 
can be shown that for a polymer with either free ends or 
forming a closed loop (without any further restrictions 
on topology) 
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where a is a microscopic length-scale (of the order of 
monomer diameter or bond length), while v w 0.59 (in 
3D) is a system independent (universal) exponent. The 
function $ includes corrections to the leading power-law 
due to irrelevant variables, which can be interpreted as 
additional length scales {N}. Keeping only the largest 
such correction for large N leads to the second part of the 
above equation. The dimensionless constants A and B 
are again system specific, while A w 0.5 in 3D is universal 
Jl5||. If the leading correction comes from, say, the argu- 
ment No/N of <E>, the corresponding length-scale grows 
as N ~ N to , with A = 1 — 1 . 

For very long polymers, the correction term in Eq. (|l|) 
is irrelevant, and the dependence of many physical prop- 
erties on the number of monomers N can be cast in terms 
of a dependence on R g |l6). For example, consider the 
force / needed to stretch a polymer between two parallel 
walls at a distance R. We can construct two dimension- 
less quantities, fRg/ksT and R/R g (T is the temper- 
ature and ks is the Boltzmann constant), which must 
be functionally related. It is thus convenient to intro- 
duce variables /' = faN u /k B T and R' = R/aN v , and 
express the force-extension relation in the form 



/' = G(R'). 



(2) 



Simple arguments |lq] can now be used to determine the 
asymptotic behaviors of G(R'): For a large stretching 



force, the distance R between the ends of the polymer 
must be proportional to N. Conversely, in a strongly 
compressed state (R <C R g ), the force must be propor- 
tional to N. These limiting behaviors can be reconciled 
with Eq. (§) only if pi 
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for R' > 1, 
for R' < 1. 
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In the large force regime this gives R ~ 
aNifa/kBT)^-^/" = aN£(N/N b ) (omitting di- 
mensionless prefactors), with Nt, = {ksT '/ fa) 1 '" . Thus, 
the polymer can be viewed as a linear sequence of 
N/Nb blobs, each of size aN^ and consisting of Nb 
monomers jl7|]. On length-scales smaller than the blob 
size, the external forces are not significant, while on the 
length-scales larger the polymer is essentially linear. 
(An analogous blob picture is also available for the 
compressed regime J_6|.) 

We employed Monte Carlo (MC) simulations to mea- 
sure such force-extension relations. Our model chains 
were composed of hard spheres of diameter 0.75a con- 
nected by "tethers" restricting the distance between ad- 
jacent spheres to be smaller than a, with no additional 
energy costs. The end-monomers were fixed to two in- 
finite parallel walls a distance R apart. An elementary 
MC step consisted of an attempt to move a randomly cho- 
sen sphere a distance 0.16a in a random direction. (./V 
such attempts constitute one MC time unit.) With such 
parameters the chain cannot cross itself, and its topol- 
ogy is preserved by the impenetrable walls. The force / 
was calculated from the probability densities of contacts 
between spheres, and the probability densities of having 
stretched tethers, as described in Ref. fll8| . 

We studied chains of lengths N = 225, 335, 500, and 
750, in both unknotted (simply connected), and knotted 
(connected via a single trefoil) states Q. Fig. [I] depicts 
two different configurations of the knotted chains, where 
the region in which the knot is "concentrated" has been 
lighter shaded. For each R and N our simulations lasted 
about 10 8 MC time units, which is considerably longer 
than the estimated Rouse relaxation time We ver- 
ified that during the relaxation period the knot is able 
to diffuse from one side to the other. Open and closed 
symbols in Fig. ^ depict the results for unknotted and 
knotted chains, respectively. 

It is important to note that the usual derivation of 
Eqs. ([!])- (|^) is not for a specified polymer topology, but 
rather for an ensemble of polymers that includes all pos- 
sible topologies. However, relations of this type are also 
likely to be valid for polymers of fixed topology, albeit 
with a different scaling function in Eq. ([!]). The collapse 
of the data for unknotted polymers (open symbols) of 
different lengths in Fig. |^ confirms this expectation. In- 
deed, the quality of the collapse indicates that subleading 
corrections are negligble for R' < 2. We have taken ad- 
vantage of this observation to construct an analytic fit 
to the scaling function G, as depicted by the solid line 
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in this figure. For R' > 2 the polymer with N = 335 
forms more than six blobs, each containing less than 55 
monomers. For such small blobs, finite size effects begin 
to appear, as we observe a roughly 10% deviation in the 
right-most data points in Fig. |[ Note that a deviation 
of this magnitude is expected from Eq. |l|, with B ~ 1. 
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FIG. 2: The scaled force /' as a function of the scaled sep- 
aration between walls R' , for chains of sizes N = 750 (cir- 
cles), N = 500 (squares), N = 335 (triangles), and N = 225 
(diamonds) monomers. Open and solid symbols correspond, 
respectively, to unknotted and knotted chains. The solid line 
depicts an analytical fit to G in this range. 

As expected, for a given /' the knotted polymer (solid 
symbols in Fig. ||) has a significantly smaller R' than its 
unknotted counterpart. However, the scaled difference 
becomes less pronounced for larger N, and the results 
for longer chains approach those of the unknotted chain. 
The absence of data collapse indicates the appearance of 
strong finite size corrections in a parameter range where 
the unknotted chains show no such effect. We would like 
to associate this feature with the emergence of a new 
size scale due to the knot. If a string with a knot is fully 
stretched, its maximal length is reduced by the size of 
the resulting tight knot. While the knots in our sim- 
ulations are far from tight, we shall still describe the 
influence of the knot as a reduction in the number of 
monomers N by the 'size of the knot' N k . If the knotted 
chain of length N is equivalent to an unknotted chain of 
length N — N k , its force-extension curves must satisfy 
fa(N - N k ) v /k B T = G[R/a(N - N k )% where G is the 
scaling function obtained before for unknotted polymers 
(solid line in Fig. ||). Naturally, this definition will not 
work with a single Nk, since our knots are not tight. How- 
ever, each solid data point in Fig. || can be moved to the 
previously obtained solid line by choosing an appropri- 
ate N k {f',N). This can be regarded as our operational 
definition of the size of the knot for a given / and N. 

If the force-extension curves for knotted polymers are 



to be consistent with the standard finite size corrections 
discussed in connection with Eq. (Q), we must have 

/' = G k (&, ^, . . . ) ~ G(R') [1 + g(R')N-^} . (4) 

Based on the numerical results, we have assumed that 
the leading scaling function G is the same for knotted 
and unknotted polymers, but allowed for different cor- 
rections, such as a new exponent For Eq. (||) to be 
consistent with our definition of knot size, we must have 
N k (f,N) = H(f')N t . The exponent t = 1 - A k gives 
the scaling of the number of monomers in a knot (in the 
absence of force) via N k ~ N*. General scaling consid- 
erations do not restrict the shape of H(f'). However, 
from the "blob picture" of a strongly stretched chain, 
we know that the polymer is essentially linear beyond 
the blob size iVj, (and consequently not knotted on such 
scale), while within a blob it is undisturbed by the exter- 
nal forces. We thus expect the knot size to be determined 
by Nf, (as if this is the entire length of the polymer), i.e. 

N k ~ N* = N t /f' t ' v , and H(f') ~ 1/ '/'*'" for /' » 1. 
Similar behavior occurs |2l| in force-extension charac- 
teristics of polymers in which a sliding constriction (or 
slip- link p2|) creates a loop, somewhat reminiscent of a 
topological constraint. In the latter, the size of the loop 
under strong tension is equal to the size of a loop in un- 
stressed polymers of size Nf,. 
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FIG. 3: The scaled knot size N k /N 0A , as a function of the 
scaled force /', for chain lengths N = 750 (circles), 500 
(squares), 335 (triangles), 225 (diamonds). 

While our range of extensions is too limited to test the 
asymptotic behavior of H(f'), we can estimate the ex- 
ponent t from the value at which the functions N k /N l 
exhibit the best collapse, i.e. are least sensitive to N. 
Fig. U depicts the optimal collapse with t = 0.4, which 
is characterized by a \ 2 value of 0.35. For t — 0.3 and 
t = 0.5 we have x — 1 5 serving as a criterion for the 
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error in t, and our estimate of t — 0.4±0.1. This is some- 
what smaller than the value that can be deduced from 
Ref. juj . The corresponding finite-size correction expo- 
nent is Afe = 0.6. We note that this is close to the best 
numerical estimate A ~ 0.56, for the dominant correction 
to scaling when considering all topologies p3j , and within 
errors, also is consistent with the estimate A ~ 0.48 ob- 
tained by field theoretic techniques [p4[. Is this more than 
simple coincidence? The standard field theory for poly- 
mers is based on an expansion around four dimensions 
that does not incorporate topological constraints. As- 
suming that the analytic continuation of this theory to 
three dimension also tells us about knots, how can such 
effects be anticipated in the perturbative expansion? If 
knots do indeed appear with sublcading sizes in 3D, their 
effects could be anticipated in corrections to scaling, in 
which case A^ = A. This conjecture is bolstered by 
A ~ A fc ~ 1 for knots confined to 2D [g, [||. 

In summary, by comparing force-extension relations of 



unknotted and knotted polymers of several lengths, we 
observe strong finite-size corrections in the latter, which 
we attribute to the knot size Nk- Scaling analysis and 
data collapse suggest a power law Nk ~ N*, with t = 
0.4 ± 0.1. Thus unlike 2D 'flat knots' §, the 3D knot 
sizes grow with the length of the polymer (although as a 
diminishing fraction of the whole length). Since a single 
knot is only weakly localized and "knows" about the size 
of the chain, it is interesting to investigate chains with 
several knots which may interact with each other. This, 
however, requires simulations with much larger chains. 
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